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Abstract 

Three-dimensional hypersonic rarefied flow about the Aeroassist Flight 
Experiment (AFE) vehicle has been studied using the direct simulation Monte 
Carlo ( DSMC ) technique. Results are presented for the transitional flow 
regime encountered between 120 and 200-km altitudes with a reentry velocity 
of 9.92 km/s. In the simulations, a five-species reacting real-gas model 
that accounts for internal energies (rotational and vibrational) is used. 
The results indicate that the transitional effects are significant even at 
an altitude of 200 km and influence the overall vehicle aerodynamics. For 
the cases considered, the aerodynamic coefficients, surface pressures, 
convective heating, and flowfield structure variations with rarefaction 
effects are presented. 

Introduction 

Over the past 20 years, increasing attention has been given to 
aeroassisted flights in which a space vehicle is decelerated and/or 


maneuvered through a planet's atmosphere by making use of aerodynamic forces 
rather than an all-propulsive system. The objective of these efforts has 
been to reduce the weight of the propellants carried by the vehicle and to 
increase the payload. Today, attention has been focused on applying this 
idea to reusable aeroassisted orbital transfer vehicles (AOTV) which would 
transfer payloads between two orbits in space. Recently, the National 
Aeronautics and Space Administration initiated the Aeroassist Flight 
Experiment (AFE) in order to provide flight data to be used in the design of 
an AOTV. In this experiment, a subscale model vehicle will be deployed from 
the Space Shuttle and accelerated by a solid rocket motor to enter the 
Earth's atmosphere down to an altitude of about 75 km. The vehicle will 
then exit the atmosphere to be retrieved by the Shuttle. Aerothermal loads 
will be constantly recorded by the AFE vehicle during this representative 
flight, which is scheduled for the mid 1990's. 1 

Today, estimation of the aerothermal loads for the design of an AOTV 
and the AFE vehicle primarily depends on numerical simulations because 
ground-based facilities fail to simulate the highly nonequilibrium 
conditions about these vehicles. After the AFE experiment, the recorded 
flight data will be used to validate these simulations. 

Analysis 

At high altitudes, the mean free path of molecules is comparable with 
the size of space vehicles, thereby invalidating the use of the Navier- 
Stokes equations because the continuum assumptions are not applicable. For 
these flows, the direct simulation techniques appear to be the most feasible 
methods from the computational point of view. 
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In this study, the direct simulation Monte Carlo (DSMC) method as 
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developed by Bird is used. This method simulates transitional flows by 
following the motion of thousands of representative molecules in the 
fiowfield. The molecular motion is computed for small time intervals. 
Intermolecular collisions appropriate to these time intervals are uncoupled 
from the molecular motion and are calculated by choosing collision partners 
on a probabilistic basis. In this study, the collisions are modeled by the 
variable-hard-sphere (VHS) model, as is the case with most contemporary 
DSMC applications. The VHS model treats the molecules as hard spheres, as 
far as the scattering is concerned, but the collision cross section depends 
on the relative speed of the colliding partners. A reacting real-gas model 
with five chemical species (0^, N^, 0, N, and NO) is used for these 
computations. Ionization is neglected on the expectation that it will not 
be significant for the present conditions. 

For computational convenience, the domain is divided into regions, and 
each region consists of a network of cells. In the three-dimensional DSMC 
applications, the cells are deformed hexahedra, and each cell is further 
divided into five tetrahedral subcells. The collision partners are selected 
from the same tetrahedral subcells (rather than from the same cells) to 
reduce the distance between the colliding molecules and to improve the 

2j 

results with regard to the conservation of angular momentum. 

Computational Results 

In this study, flow about the AFE vehicle has been investigated. The 
shape of the AFE vehicle is that proposed by the NASA Johnson Space Center 
and is sketched in figure 1. The vehicle basically consists of a forebody 
aerobrake, a hexagonal experimental carrier, and a solid rocket motor which 
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is ejected prior to atmospheric entry. The axis of the solid rocket motor 
is inclined at a 1 7° angle with respect to the incoming freestream. 

The angle of incidence remains constant for the present solutions. The 
wind vector was assumed to be parallel to the x-axis (fig. 1), which 
corresponds to zero incidence for the present coordinate system. (This is 
equivalent to an angle of incidence of 17° for the AFE body axes system in 
which the x-coordinate is coincident with the solid rocket motor axis.) 

The forebody of the AFE vehicle is generated by a 60° elliptic cone 
raked at a 73° angle. Its nose is blunted in the shape of an ellipsoid, and 
the aft skirt is appended such that at each meridional plane it is tangent 
to the cone by a circle of fixed radius. The base length of the vehicle is 

5 

4.25 m. Cheatwood et al. describe mathematically the shape of the AFE 
forebody and have written a computer program to generate the forebody 
surface coordinates. Their program has been Incorporated in these 
computations to describe the forebody surface grid. The coordinate system 
used in this study is located at the tip of the elliptic cone such that the 
x~axis becomes parallel to the freestream velocity, and the x-y plane is 
located on the symmetry plane, as shown in figure 1 . According to this 
system, the geometric stagnation point of the AFE vehicle lies at x = 0.3 m» 
y = 0, z = 0. 

In general, calculations for high altitudes need to be performed on 
large computational domains where the cells can be quite large because the 
molecular mean free path, and consequently the disturbance field in front of 
the body, is large. This fact is reflected in the present calculations. 

For example, the computational domain in front of the body starts at 
x = -4.7 m for the 1 20-km-altitude case, at x = -10.3 m for the 140-km- 
altitude case, and at x = -18 m for the 150-km and higher altitude cases. 



The 1 30-km-alt itude computations show that the afterbody portion does not 
affect the vehicle aerodynamics, hence for the 140-km and higher altitudes, 
only the forebody domain is used. For all the calculations, only half of 
the geometry is considered because of the plane of symmetry. 

Figure 2 shows the computational grid used for the 1 20-km-altitude 
case. In this figure, both the cells and subcells are shown on the outer 
freestream boundary. However, on the plane of symmetry, only the cell 
structure is drawn for clarity. For all cases considered, the cells and 
subcells are exactly matched with the neighboring ones throughout the 
computational domain. For the computations where both the forebody and 
afterbody are considered (i.e., the 120 and 1 30-km-altitude cases), a total 
of 21 regions with 1 ,820 cells is used, and these flows are simulated using 
about 55,000 molecules. The forebody domain is divided into 6 regions with 
600 cells, and the 140-km and higher altitude cases are simulated using 
about 20,000 molecules. In all cases, the AFE forebody surface consists of 
100 triangles, or in other words, 50 distorted rectangles with 10 divisions 
in the longitudinal direction (as can be seen from fig. 2) and 5 divisions 
in the spanwise direction. 

At the freestream boundaries, atmospheric conditions appropriate for 
the particular altitude are used. Table 1 lists the atmospheric conditions 
for altitudes between 120 and 200 km. The uniform forebody temperature i3 
based on free-molecule radiative equilibrium heat transfer to the stagnation 
point. These values are tabulated in Table 2. The AFE vehicle is insulated 
with very low conductivity materials. The afterbody wall temperature is, 
therefore, assumed to be at a uniform 300 K for all cases. At the body 
surface, full thermal accommodation with diffuse reflection for the gas- 
surface interaction is assumed. 
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Once steady state was achieved in the simulations, computations were 
continued for an additional 8,000 time steps, and samples were taken every 
other time step to arrive at the time-averaged results. Examples of the 
time-averaged results presented in figures 3~12 are based on these 4,000 
samples. In figures 3~6, surface pressure and heat transfer rate 
distributions are presented. Both pressure and surface heating have their 
maxima slightly below the geometric stagnation point, are fairly uniformly 
distributed along the forebody surface (this was one of the design aims), 
and quickly decrease along the skirt. 

Figures 7~9 show various flowfield contours in the x-y plane for the 
120-km case. In these figures, the vehicle is half a cell size larger in 
all directions. The distortion in the picture is greatest in the afterbody 
where rather large cells are used. Clearly, the flowfield density reaches 
its maximum at the body surface, whereas the maximum translational (kinetic) 
temperature occurs somewhere away from the surface in the upstream 
direction. On the other hand, the internal temperature shows a rather 
complicated picture and is discussed in the next section. 

Figures 10-12 and Table 3 present the aerodynamic coefficients 
(C = drag coefficient, C - lift coefficient, L/D = lift-to-drag ratio) at 

D Li 

selected altitudes for an angle of incidence of 0° (using the present 

coordinate system shown in figure 1). The forces are normalized with 

2 

respect to 1/2 p U A _ where p and U are the freestream density and 

r CD 00 p 0 J r 00 00 J 

2 

velocity, respectively, and A f,, the reference area, equals 14.1 m . These 
figures and the table also show the calculated free-molecule and modified 
Newtonian results, along with experimental wind-tunnel data. The 
experiments were conducted in the NASA Langley Research Center Mach 10 air 


6 


and Mach 6 CF^ (freon) wind tunnels using high-fidelity models. ^ The wind- 
tunnel values presented in Table 3 and figures 10-12 are recommended values 
for zero incidence based on the Mach 10 air and Mach 6 CF^ results (personal 
communication from Chris Cerimele, NASA Johnson Space Center). Clearly, the 
DSMC results approach the free-molecule limit very slowly at higher 
altitudes, and even at an altitude of 200 km, the flow is not completely 
collisionless. Prior to this study, it was generally acknowledged that 
free-molecule flow existed for the AFE vehicle for altitudes near 150 km, 
but this study shows that the transitional effects are significant at these 
altitudes and influence the overall aerodynamic coefficients. These results 
have important implications for the interpretation of aerodynamic 
coefficients extracted from flight measurements under rarefied conditions. 
Figures 10-12 also contain the results of the Lockheed bridging formula 
which empirically connects the axial and normal aerodynamic force 
coefficients between the continuum and free-molecule limits. This is 
accomplished with a sine-square function by assuming continuum flow at a 
Knudsen number Kn =0.01 and free-molecule flow at Kn =10, which 
correspond to altitudes of 90 and 150 km, respectively. The bridging 
formula results are plotted to show the general trend even though its 
results are erroneous for the conditions considered in the present study. 
Figures 10 and 11 also show effects of the pressure and shear forces on the 
total drag and lift coeff icients. Clearly, the pressure contributions are 
fairly constant at all altitudes and are approximately equal to the measured 
continuum wind-tunnel values. However, the shear contribution increases 
with altitude, and in effect, increases the drag and reduces the lift. The 
pressure coefficient at the stagnation point is 1.944 for the 1 20-km case 
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and 2.033 for the 200-km case. The corresponding limiting value for the 

continuum regime is 1.839, and the free-molecule result is 2.051. 

One of the design aims of the AFE was to have the vehicle trim at IT 0 

angle of attack. At this attitude, therefore, the pitching moment should be 

zero. The computed pitching moment coefficients with respect to the center 

-2 

of gravity are less than 10 , and unlike the drag and lift, are expected 

to be sensitive to the surface grid chosen. Considering the relative 
coarseness of the grid and the scatter due to the small sample, it is more 
meaningful to report the equation of the line along which the resultant 
force acts. For the 1 20-km case, this line is y = 0.07x - 1.031 according 
to the coordinate system shown in figure 1 and varies smoothly with altitude 
to y = 0.01 3x - 0.984 for the 200-km case. 

Discussion 

The computations in this study have been performed using a general 3-D 

DSMC computer program written by G. A. Bird. This program has already been 
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applied to various wedge flow problems by the present authors. Three- 
dimensional DSMC applications have recently become possible through 
advancements in computer technology, and various cylinder/cone problems have 

g 

been studied by other groups as well. However, this is the first time that 
a general 3~D DSMC code has been used to describe the flow about a full- 
scale space vehicle. 

The first problem investigated in this study was the case where the AFE 

vehicle is represented by an idealized axisymmetric geometry. This case was 

previously studied using an axisymmetric DSMC code for an altitude of 1 30 km 
q 

and below. Therefore, initial application of the 3~D program concentrated 
on the axisymmetric flow problem at 1 30 km because higher altitude flows 
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are, in general, easier to compute. The calculated results were in very 
good agreement with those of the axi symmetric code computations. The grid 
used in the 3 -D computations (8 divisions along the surface and 16 divisions 
in the normal direction) was coarser than the one used in the axisymmetric 
program (20 divisions along the surface and MO divisions in the normal 

direction), and the cell structuring was different. Agreement of these 

\ 

results, consequently, is an indication of grid independency for the 
computations. For this case, the extent of the computational domain was 
also checked. When the domain was started 5 m in front of the body (where 
density ratio p/p^ = 1.5 and translational temperature T = 17,000 K) 
instead of 15 m (where p/p = 1.0 and T = 1,500 K) , the results showed 

03 UP 

almost no difference. 

Subsequently, the 3~D program was applied to the three-dimensional flow 
about the actual AFE vehicle at 130 km, first without and then with the 
afterbody. It should be noted that the three-dimensional shape of the AFE 
vehicle is exactly generated without any approximations. In the 
computations, only the forebody and the experimental carrier were included 
because the solid rocket motor is ejected during entry near 1 30 km. Above 
this altitude the afterbody has no effect on the aerodynamics. The 
experience accumulated on the 130-km altitude case has been used for the 
other altitude computations, and only one computation has been performed for 
each of those cases. Nevertheless, in all of those computations, the basic 
DSMC requirements are satisfied. Namely, the cell sizes are smaller than 
the local mean free path, the number of molecules in each cell is of order 
10, and the time step used to move the molecules is smaller than the mean 
collision time per molecule. Past experience indicates that these results 
will not differ if finer grids are used. 
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In analyzing these results, one has to remember that these flows are in 
a highly nonequilibrium state and that these macroscopic properties are 
obtained by averaging over various classes of molecules whose distribution 
functions are widely separated in velocity space. Therefore, it is not 
possible to properly interpret these kinetic temperature contours entirely 
from the viewpoint of classical thermodynamics. 

The complicated internal energy distribution shown in figure 9 is due 
to the three distinct classes of molecules present in the flowfield. First 
of all, at these altitudes there are quite a few freestream molecules with 
low internal energies crossing the entire domain without any collisions. 
Second, there is the effect of intermolecular collisions. Although there 
are very few collisions, the internal energies of the molecules which have 
undergone a collision are altered drastically because of the energy exchange 
with the high translational energies of the gas molecules. Third, there 
are molecules that have been reflected from the body surface. These 
molecules have low internal energies because the surface temperature is low, 
and it is assumed that these molecules are reflected with full thermal 
accommodation. The relative abundance of these three classes of molecules 
varies with location in the flowfield, and their mixing produces a highly 
nonequilibrium flow. Figure 9 shows the effects of this on the macroscopic 
flowfield properties. As the flow approaches the body along the stagnation 
streamline, the average internal energy first increases as a result of 
intermolecular collisions and then decreases because of the increase in the 
concentration of sur face-reflected molecules. As the molecules move around 
the body, they have more opportunity to have collisions, and since the 
relative abundance of the collided molecules increases, so does the average 
internal energy. In the wake region, the molecular number density is 
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extremely small due to the "shadow effect" of the forebody. Consequently, 
this region has basically one class of molecules, and they are the ones 
which have undergone a collision. Because this class of molecules has high 
internal energies, the average internal energies in the wake region tend to 
be higher. 

As the altitude increases to 200 km, the internal temperature picture 
changes entirely. This is not surprising, because at this altitude the 
freestream temperature is higher, the body surface temperature is cooler 
(see Tables 1 and 2), and more importantly, the number of intermolecular 
collisions is much lower. In other words, the flowfield has basically two 
classes of molecules, namely, freestream and surface-reflected molecules. 
With full thermal accommodation, the surface-reflected molecules have low 
internal energies due to the cool surface temperature, and naturally, their 
concentration is high near the body. Consequently, the average internal 
energy in the forebody region continuously decreases as the flow approaches 
the body. 

For 10 km/s flows under highly rarefied conditions, the assumption of 
full thermal accommodation is almost certainly incorrect. As stated before, 
for most cases considered in this study, there are quite a few freestream 
molecules which hit the body surface with very large translational energies. 
According to the full thermal accommodation assumption used in this study, 
all of the energy of the molecules is converted to surface heating, and the 
molecules are reflected with low internal energy that is characteristic of 
the surface temperature. However, this assumption is probably unrealistic, 
and it is quite likely that some of the translational energy appears in 
excited internal modes during surface reflection. Actually, this phenomenon 
has already been observed on the Shuttle Orbiter. 1 ^ Unfortunately, the 
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mechanism underlying this complicated interaction and the nature of the 


energy exchange are not known. These questions need to be clarified in 
future experiments. Until then, the internal temperature results of this 
study should be interpreted with reservation. At lower altitudes, this 
phenomenon is less significant because the molecules v/hich hit the surface 
have less translational energy due to intermoleeular collisions. Even 
though the nature of the gas-surface interaction could significantly alter 
the internal temperatures, it should not affect the translational 
temperatures and the aerodynamic coefficients because the translational 
energy is most likely fully accommodated to the surface. 

Summary and Conclusions 

This study applies the direct simulation Monte Carlo technique to 
assess the three-dimensional flow about the Aeroassist Flight Experiment 
(AFE) vehicle. For altitudes of 120 to 200 km, results show that: 

(1) The flow approaches the free-molecule limit very gradually at 
higher altitudes, and even at 200 km, the flow is not completely 
collisionless . 

(2) Maximum surface pressure and heating occur slightly below the 
geometric stagnation point of the vehicle. 

(3) At altitudes of 1 30 km and above, the effect of the AFE afterbody 
on the vehicle aerodynamics is negligible. 

(4) Dissociation is negligible. 

This is the first time that the actual three-dimensional flow about the 
AFE vehicle has been studied at altitudes of 120 km and above. For the 
cases investigated, radiation is negligible and is not considered. At these 
altitudes, the aerothermal loads on the vehicle are small, and the peak 
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loads and heating occur at much lower altitudes. This should not overshadow 
the significance of the present study because these results are the first 
3-D solutions to the AFE flow problem in the transitional regime, and they 
show that the transitional effects are significant even up to an altitude of 
200 km. Results such as these are essential for the proper interpretation 
of the high-altitude aerodynamic measurements. 
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Table 1. Freestream conditions. 


Mole fractions 


Altitude, 

km 

Density 

kg/nr 

x10 9 

, Temper atur 

K 

x o 2 

X N. 

x o 

Mean free path, 
m 

200 

0.3 

1026 

0.031 

0.455 

0.514 

197.00 

170 

0.9 

892 

0.044 

0.548 

0.408 

75.61 

150 

2.1 

733 

0.055 

0.616 

0.330 

30.93 

140 

3.9 

625 

0.062 

0.652 

0.286 

16.92 

1 30 

8.2 

500 

0.071 

0.691 

0.238 

7.72 

120 

22.6 

368 

0.085 

0.733 

0.183 

2.69 

Table 

2. AFE 

forebody wall 

temperatures 

used in 

the computations. 


Altitude, 

km 

T 

wall, 

K 

200 

230 

170 

300 

150 

370 

140 

425 

130 

500 

120 

650 
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Table 3. Aerodynamic coefficients. 



Altitude , 
km 

C D 

C L 

L/D 

Calculated 




(Free molecule) 

1 .97 

0.019 

0.010 

200 

1 .95 

0.025 

0.013 

170 

1 .93 

0.035 

0.018 

150 

1 .90 

0.049 

0.025 

1 40 

1.87 

0.064 

0.034 

1 30 

1 .84 

0.086 

0.046 

120 

1 .80 

0.125 

0.070 

(Modified Newtonian) 

1.24 

0.338 

0.274 

Experimental 




(Wind tunnel) 

1.35 

0.381 

0.281 
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Experimental carrier 



Figure 1. The AFE vehicle. 
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Freestream 

boundary 



Plane of 
symmetry 


Figure 2. Computational grid 
(Altitude = 120 km). 
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Figure 3. Surface pressure distributions. 
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(b) Spanwise direction 

( s 2 is measured along the surface from the stagnation point). 


0.6 



Maximum 
h 2.22 N/m 2 


Figure 4. Surface pressure contours 
(Altitude = 120 km). 
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Figure 5. Surface heat transfer rate distributions. 






Figure 6. Surface heat transfer rate contours 
(Altitude = 120 km). 
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Figure 7. Flowfield density contours at the symmetry plane 
(Altitude = 120 km). 
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Figure 8. Flowfield translational temperature contours at the symmetry plane 
(Altitude = 120 km). 
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Figure 9. Flowfield internal temperature co 
(Altitude = 120 km). 
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Total drag (DSMC) 
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Figure 10. Drag coefficient variation with altitude (FM denotes free - molecule 
value; the wind-tunnel value is a composite of Mach 10 Air 
and Mach 6 CF 4 data). 


Wind-tunnel Pressure contribution 
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Figure 11. Lift coefficient variation with altitude (FM denotes free - molecule 
value; the wind-tunnel value is a composite of Mach 10 Air 
and Mach 6 CF 4 data). 
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Figure 12. Lift /Drag variation with altitude (FM denotes free - molecule 
value; the wind-tunnel value is a composite of Mach 10 Air 
and Mach 6 CF4 data). 
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